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Linear stability of an isothermal, pressure-bounded, self-gravitating gas slab which is 
gravitationally coupled with the background weakly interacting massive particles (WIMPs) 
O |. is investigated. Analytic dispersion relations describing such a configuration are derived. 

Two novel, distinct oscillatory modes are found. Astrophysical implications of the results are 
discussed. 

Subject headings: dark matter — gravitation — instabilities — ISM: clouds — ISM: kinematics 
and dynamics — stars: formation 
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It is widely accepted that the formation of stars occurs through gravitational collapse of the interstellar 
gas clouds. One of the possible mechanisms of triggering the star and/or stellar cluster formation process is 
a high-velocity (supersonic) cloud-cloud collisions. In a such collision a dense gaseous slab is formed with 
two plane-parallel shock fronts propagating away from the interface (Usami et al. 1995). Then the slab 
grows in mass becoming unstable against gravitational instability which causes its fragmentation. Newly 
formed the gaseous clumps, in turn, collapse further and evolve into stars and/or stellar clusters. 
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Introduction 



The structure of a self-gravitating, isothermal shock pressure-bounded slab of gas which usually is 
formed in the cloud-cloud collisions was investigated by Ledoux (1951) and later by Elmegreen & Elmegreen 
(1978) among others. Lubow & Pringle (1993), LP93 thereafter, have also derived an analytic dispersion 
relation for the linear perturbations and investigated, in detail, the full unstable branch in such a slab of 
gas. 

Possible existence of WIMP matter, one of the form of dark matter which itself is a dominant mass 
component of the universe, is strongly motivated both by standard models of particle physics and cosmology. 
It has been established that the mass density associated with the luminous matter (stars, hydrogen clouds, 
x-ray gas in clusters, etc.) cannot account for the observed dynamics on galactic scales and above (Trimble, 
1987), therefore, revealing the existence of large amounts of dark matter in this universe. It has been 
argued that the dark matter could be anything from novel weakly interacting elementary particles (massive 
neutrinos, neutralinos) to baryonic matter in some invisible form (brown dwarfs, primordial black holes, 
cold molecular gas, etc.). 

Recently, it has been shown (Tsiklauri 1998) that classical Jeans theory of the stability of self-gravitating 
interstellar gas clouds is significantly modified by taking into account the presence of background weakly 
interacting massive particles. It has been found (Tsiklauri 1998) that the presence of WIMP matter yields 
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an unavoidable reduction of the Jeans length, the Jeans mass and the Jeans time (time-scale of the collapse 
via gravitational instability). 

Resuming aforesaid, and in the light of recent results of Tsiklauri (1998), a revision of a conventional 
theory (e.g. LP93) of the stability of self-gravitating, isothermal shock pressure-bounded slab of gas against 
linear perturbations in the case of presence of WIMP matter, or generally speaking any type of microscopic 
non-baryonic dark matter that is coupled with the slab only via gravitational interaction, seems to be of a 
considerable importance for full understanding of the star and/or stellar cluster formation process. 



2. The model 

Equations that govern dynamics of two self-gravitating fluids (a slab of gas and WIMP matter) 
inter-coupled only via gravitational interaction can be written in the following way: 
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where notation is standard: Vi, Pi, pi and <j> denote velocity, pressure, mass density and gravitational 
potential of the fluids. The subscript i = G,W denotes two components, baryonic gas and WIMP matter 
respectively. Now, writing every physical quantity, for brevity commonly denoted by f(x,z), in a form of 
f(x, z) = f" + f'(z) exp[i(wi — kx)] and doing usual linearization of the Eqs.(l)-(4) we can obtain following 
a set of coupled differential equations for the perturbations 
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Here, C s g and C s w denote speeds of sound for the baryonic gas and WIMP matter respectively. It is 
worthwhile to note that Eqs.(5) and (6) resemble closely to the equations describing coupled pendulums 
which arise in the non-modal study of linear perturbations in shear flows (cf. Chagelishvili, Rogava & 
Tsiklauri 1996, 1997). The fundamental (normal) oscillatory modes of the Eqs.(5) and (6) are given by the 
following expression: 
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are the eigenfrequencies of the system. Also, a G = Ai:Gp G /C 2 G and aw = ^Gp w /C 2 w . 



A general solution of the Eqs.(5)-(7) can be readily found. Following LP93 we consider only symmetric 
modes dP' G /dz = dP w /dz = at z = 0. There are two distinct, general solutions of the Eqs.(5)-(7) 
depending on the sign of q 2 _ . 

For the q 2 _ > we obtain 

P' G = A- cos((7_z) + A + cos(q + z), (11) 

2 2 2 2 

Pw = q -~ q ° A- cos(g_z) + q+ ~ qa A+ cos(q+z), (12) 
aw aw 

$ = Ccosh(fcz) - (p G + q -~ q2G f3w) ]2 A ~ 2 cos(g_*)- 
\ a w ) k 2 + q± 



Pg + + U Pw tt-^ cos((Z + z), (13) 
aw J k + q+ 

here we introduced notation Pg = ^G/C 2 G and Pw = 4:TrG/C 2 w . 
For the q 2 _ < we obtain 

P' G = A- cosh(y / -g 2 _ z ) + A + cos(q + z), (14) 

2 2 2 2 

P'w = ~ q ~~ qG A- coshU Pq^z) + q+ ~ qG A+ cos(q + z), (15) 
aw v aw 

cj) 1 = Ccosh(fcz) - (p G + ~ q ~~ qG pw) jy— — g- cosh(,/^"z)- 



v avi/ / K — qz. 

(p G + ^^Pw) ™ S (q+z), (16) 
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Note that in the case of presence of the background WIMP matter there are two fundamental (normal) 
modes q± (see Eq.(8)). Whereas, in the case with no WIMP matter there is only one oscillatory mode, 
namely qa (see Eq.(32) in LP93). This is the presence of WIMP matter that causes appearance of a novel 
unstable mode described by Eqs.(14)-(16). 

As our next step towards derivation of the dispersion relation for linear perturbations, following LP93, 
we impose relevant boundary conditions. The first boundary condition is that the Lagrangian pressure 
perturbation is zero at the slab boundary, namely that 

dP° 

Pi + ^=0, (17) 

at z — a (with a being slab half-thickness), where £ vector is the Lagrangian displacement vector of the 
perturbation and i = W,G. In this case, the displacement is simply related to V' by £ = —iV'/u. Using 
relation VPj° = — p° V<^>° = —ATtGap®{p G + p w ) for the unperturbed physical quantities, Eq. (17), at z = a, 
reduces to the two following conditions 

^=^G(p G +p° w )a, (18) 
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The other boundary condition is the one for <j>' (LP93): 

#L dxjj 



dz dz 



= 47rG(p° G +p° w %, 



evaluated at z = a, where <j>' e is the perturbed gravitational potential in the region \z\ > a and is given by 
4>' e = (f>'(a) exp[— k(\z\ — a)], where it has been assumed that k > 0. Thus, we obtain 



dz 
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2.1. The case when q 2 _ > 



Applying boundary condistions (18) and (19) to the general solutions in the case when q 2 _ > 
[Eqs.(ll)-(13)] yields 
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Another condition which relates A + and C can be obtained using z-componcnt of the Eq.(2) for the 
perturbations in baryonic gas and £ = —iV'/u) (LP93). Doing this yields 
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Now, combining Eqs.(20) and (21) allows us to obtain the dispersion relation for the perturbations for the 
case when q 2 _ > 
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where b^ = 1 — D + p G ak, b 1 = — x(l — D-p G ak), b 2 = D + p G q + a, b 2 = — xD_p G q_a, 
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2.2. The case when q 2 _ < 



We now apply boundary conditions (18) and (19) to the general solutions in the case when q 2 _ < 
[Eqs.(14)-(16)]. Thus we obtain 
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Here we have introduced following notation: 
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Again using the condition which relates A+ and C that is obtained using z-componcnt of the Eq.(2) 
for the perturbations in baryonic gas and £ = —iV'/oj, we get 
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Finally, using Eqs.(23) and (24) we to obtain the dispersion relation for the perturbations in the case when 
<Z 2 <0 
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3. Discussion 

In this paper we have derived two dispersions relations for linear perturbations in a pressure-bounded, 
self-gravitating gas slab which is gravitationally coupled with the background weakly interacting massive 
particles. In the original paper LP93 authors found that the slab alone is unstable at sufficiently long 
wavelength and the growth rate is of the order of \/Gp, with p being mass density of the slab, but at 
high external pressures the nature of the instability is quite different from the classical Jeans instability. 



- 6 - 



Moreover, they have developed an analytic model which reproduces numerical results of Elmegreen & 
Elmegreen (1978). Interesting result of these authors is that the critical wavenumber for the onset of the 
instability is always of the order of the slab thickness, regardless of the level of self-gravity and the external 
pressure. 

Present analysis revealed that taking into account presence of the putative WIMP dark matter 
significantly modifies results of LP93 with major changes being appearance of two distinct fundamental 
(normal) oscillatory modes q± and the existence of an additional unstable mode when < 0. Obtained 
dispersion relations Eqs.(22) and (25) need detailed investigation in various limiting cases. Especially 
interesting would be doing analysis of the low-frequency modes, for which LP93 found an unexpected result 
that the critical wavenumber for the onset of the instability is always of the order of the slab thickness. 
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